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ABSTRACT 



We describe a one-parameter family of models of stable spherical stel- 
lar systems in which the phase-space distribution function depends 
only on energy. The models have similar density profiles in their 
outer parts {p oc r~^) and central power-law density cusps, p oc r^"'', 
< ?7 < 3. The family contains the Jaffe (1983) and Hernquist 
(1990) models as special cases. We evaluate the surface brightness 
profile, the line-of-sight velocity dispersion profile, and the distribu- 
tion function, and discuss analogs of King's core-fitting formula for 
determining mass-to-light ratio. We also generalize the models to a 
two-parameter family, in which the galaxy contains a central black 
hole; the second parameter is the mass of the black hole. Our models 
can be used to estimate the detectability of central black holes and 
the velocity-dispersion profiles of galaxies that contain central cusps, 
with or without a central black hole. 



1 Introduction 

Despite the increasing sophistication of numerical models of stellar systems, analytic or 
semi-analytic models still have many uses: they offer insight into the behavior of more 
realistic and complicated models, they provide tests of and initial conditions for iV-body 
experiments, they are free of noise and spurious numerical relaxation, and in some cases 
they fit observed stellar systems remarkably well. 

An additional motivation for the models described below is recent high-resolution 
ground-based and Hubble Space Telescope photometry of elliptical galaxies and bulges. 
These observations typically fail to resolve the central regions into constant density cores, 
even at a resolution of 0'!l. Specific examples include M32 and M87 (Lauer et al. 1992a, b). 
HST observations of a sample of more than 30 elliptical and SO galaxies by our group (Lauer 
et al. 1993) establish that, at a resolution of < 0'!l, the surface brightness profile I{R) is 
usually best approximated by R '', with < 7 < 1. Most previous dynamical studies of 
galaxies have focused on models having a constant density core, that is, an analytic surface 
brightness that can be expanded as a Taylor series in the form I{R) oc 1 — AR^ -|- . . .. Thus 
the HST observations demand new dynamical models to predict the kinematic properties 
of the central regions and to assess whether massive central black holes are implied by the 
observations. A natural first step is to seek appropriate analytic models. 

Among the simplest stellar systems are isotropic spherical systems, in which the grav- 
itational potential $ and the density p are spherically symmetric and the phase-space dis- 
tribution function depends only on the energy per unit mass, / = f{E) where E = -|-$ 
(e.g. Binney and Tremaine 1987, hereafter BT). We shall focus on the case in which the 
mass-to-light ratio T is independent of radius, so that the emissivity is p(r)/T. 
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Two of the most successful analytic models for elliptical galaxies and the bulges of 
disk galaxies are the Jaffe (1983) model, 

and the Hernquist (1990) model, 

'^'^ = 2.r{l + r)3 • (2) 

Both models are normalized so that the total mass M = Att J p(r)r'^dr is unity. The 
Jaffe and Hernquist models have similar profiles in their outer parts (p oc r~^); both also 
have central density cusps, but they differ in the strength of the cusp (p oc and 
respectively). 

In this paper we describe a one-parameter family of isotropic spherical stellar systems, 
which we call ry-models, that preserves most of the analytic simplicity of the Jaffe and 
Hernquist models, and which includes them as special cases. The parameter rj describes 
the strength of the central cusp, which diverges as /? oc r^~^ for < r] < 3. We describe 
these models in §2. In §3 we generalize the models to a two-parameter family containing 
a non-luminous central mass ("black hole" for short); the second parameter is the mass of 
the black hole. Sections 4 and 5 contain discussion. 

The work described in this paper has been anticipated to varying extents by several 
authors. Hernquist (1990, eq. 43) briefly describes a three-parameter family that includes 
the ry-models, but does not discuss the properties of any systems except the Jaffe and 
Hernquist models. Saha (1993) writes down the potential and density distributions of 
the ry-models, but focuses on their use as basis functions for solving Poisson's equation 
rather than as models of stellar systems. As this paper neared completion, we received a 
preprint from Dehnen (1993) that independently derives many of the properties of ry-models 
discussed in 52, but does not discuss models with a central black hole. 



2 Models 

We choose units in which the total mass M of the stars and the gravitational constant G 
are both unity. 

We deflne a family of spherical stellar systems called "ry-models" by the density dis- 
tribution 

Pr,(r) = ^- — -— , 0<?7<3. (3) 

We must have r\ > we show below that models with > 3 (density approaching zero 
at the origin) are unphysical. All r\ models have p oc as r — )• oo. Jaffe's model 
(eq. 1) corresponds io i] = \ and Hernquist's model (eq. 2) corresponds to = 2. The 
density diverges near the origin for all models with < 3; only models with i] = Z have 
constant-density cores. 

Saha (1993) points out that the model with 77 = has a potential corresponding to a unit point 
mass (see eq. 5 below). Thus it could represent a stellar system with luminosity density oc 1/[»'*(1 -|- »')], 
negligible mass-to-light ratio, and a central black hole. 
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Mass distribution 



The mass interior to radius r is 



r r ■ 

M^(r) = 47r pr,(r)dr = ^ 



(4) 



Potential The gravitational potential is 
^ , , f°° M^{r)dr 1 



where we assume $(r) 
finite for ?7 > 1, $,,(0) 



1] — 1 



(1 + r)*?-! 



ln(l + 1/r) 



T] = 1, 



(5) 



as r — > oo. The depth of the potential well near the origin is 
— {rj — I)""', while for < r] < 1 the potential near the origin 



diverges, — )• —(1 — rj)' 



Velocity dispersion We assume that the phase-space distribution function (df) 

depends only on energy, which implies that the velocity-dispersion tensor is isotropic. 

— 1/2 

Thus the radial velocity dispersion v"^ (r) satisfies the hydrostatic equilibrium equation, 

dr dr ' 

with boundary condition pv^ = as r — )• oo. Since p and d^/dr are non-negative, v^(r) is 
always positive. Thus we have 



V 



2 ^ 



Mr,{r)pr,{r) 



dr 



.3-.(i + ( ^ + _1 ^ + !_ 

^ ^ \2r]-4: 2r]-3 r]-l 2r] - 1 2r] 

_|_ ^^5-»? 4r''(l -|- r)^"'' 3r'^~^^(l -\- rY~'^ 

2r] -4 ^ 2r]-3 r] - 1 



+ 



2r] -1 2r] 



(7) 



This result is not valid for the four cases = |, 1, |, 2: 

^^^(r) =1 - 2r - 9r^ - 6r^ 6r^(l rf ln(l 1/r), 

<3/2(0 =r'/'(l + r)"/' - r^/'(l + r f' [41n(l + 1/r) - f ] 

- 6r5/2(l + rfl^ + 2r^/2(l + r)^/^ - \r'l\l + r)-l/^ 

^^(r) =r(l + rf [ln(l + 1/r) - ff] +4r2(l + r)2 - 3r^(l + r) + |r^ - ir^(l + r)"^ 

(8) 
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These formulae can be difficult to evaluate at large r because of near-cancellations 
between the terms. An alternate expression is the asymptotic series 



,,{r) = r'-^{l + r) 



oo 



(_2^-l)...(-2,y-j) 



1 



J! 



(5 + j)r5+i 



(9) 



Thus as r 
As r - 



oo, V. 



(r) 



2(r) 



~^ for all rj. 



r(- 



"7(4 
- In r ■ 



2^), 



12 



,3 — 77 



2j7-4 2j7 



— + — 
3 ^ 



2j?-1 



+ 



2»? 



for < ?7 < 2, 
for ?7 = 2, 

for 2 < ?7 < 3. 



(10) 



Notice the interesting range of behavior: the mean-square radial velocity diverges as 
when r] is near zero, then diverges less and less rapidly as rj increases towards 1. Once 
?7 > 1, the mean-square velocity converges to zero, converging more and more rapidly until 
?7 = 2, when it converges as r. Above rj = 2 the mean-square velocity converges to zero 
less and less rapidly, until at = 3 it is asymptotically constant. Thus the central velocity 
dispersion is constant and non-zero for two distinct models: 



!,,(r = 0) 



30 



for r] 
for r] 



1, 
3. 



(11) 



The inner parts of the model with rj = 1 resemble the singular isothermal sphere, while 
the inner parts of the model with rj = 3 resemble the central core of the non-singular 
isothermal sphere or other models with constant-density cores such as King models (King 
1966; BT). 

The qualitative reasons for the behavior of the central dispersion are worth describing. 
For T] < 2 the dispersion at a given radius is comparable to the circular speed Vc at 
that radius, given by vl = M{r)/r r''"^ which is the natural result expected from 
dimensional analysis of the hydrostatic equilibrium equation (replace d/ dr by — 1/r in eq. 
6). The behavior of models with t] > 2 is more subtle, because the mean-square velocity 
and density are dominated by different parts of the energy distribution: the density is 
dominated by low-energy stars with mean-square velocity comparable to the square of the 
circular speed, v"^ oc r^~^ ^ while the velocity dispersion is dominated by high-energy stars 
with velocities of order unity. Thus the pressure pv"^ is constant and order unity as r — )• 0; 
since p oc r^~^ the mean-square velocity v"^ oc r^~^ . 

ri/2 



The radial velocity dispersion v"^ is plotted in Figure 1 for several values of rj. 



Kinetic and potential energy The total kinetic and potential energies are 

/>oo 

Trj = Qtt / PT]{r)v^ Jr)dr 
Jo ' 

BOO 

Wr, = 27r / pr,{r)^r,{r)r^ dr. 
Jo 



(12) 
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— 1/2 . . 

1. Radial velocity dispersion function of radius (eqs. 7 and 8), for models with rj = 0.1 

(dash-dot curve), ^, |-, |- (unlabelled dashed curves), and 1, 2, 3 (solid curves). All models have 

(5r) "'/^ as r ^ oo. 

Substituting the hydrostatic equihbrium equation in the form pv^ = Mpdr/r"^ into 
the first of these equations, and the equation for the gravitational potential, $ = 
— M dr/r"^ ^ into the second, it is straightforward to confirm the virial theorem, 

BOO 

Wrj = -2Trj = -Att prj{r)Mrj{r)r dr. (13) 
Jo 
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Evaluating the integral for ry-models we have 



1 

2(2^- 1) 



Wr, = -2T,, = - ^^^^ foTf]>l; (14) 



for 1] < ^ the energies are infinite. 



Distribution function Jeans' theorem ensures that any DF that depends only on 
energy is a steady-state solution of the coUisionless Boltzmann equation. The density is 
determined by the equation 

p = 2'/\ r f{e){4,-ey/'de, (15) 
Jo 

where the DF / > is the mass per unit volume of phase space, e = — -E, and tp = — $ 
(the quantities e and tp are defined so that we work with non-negative variables: tp^r) is 
positive for all ry-models at all radii, and e is positive for all bound stars). 

Equation (15) yields a simple proof that models with rj > 3 are unphysical: since 

/ > 0, dp/dtp = 2^/^7r /J^ /(e)(V' ~ e)~^^'^de is positive; since dtp / dr = —GM{r) I r'^ is 
negative, dp/dr = (dp / dtp) / (dtp / dr) must be negative. The ry-models have p oc r^~^ near 
the origin, so models with rj > 3 are unphysical if the DF depends only on energy. 

The inverse of equation (15) is Eddington's formula (BT) 

r( . _ 1 d r dp dtp 

•'^'^ 2^/^7r^deJo dtP(e-tpy/^' ^ ' 

Integrating this equation by parts yields a simpler form: 

1 dj^ p dtp 

•^(^) = ^^7^y^ ^e-tpyi^' ^^^^ 

where we have used the fact that dp /dtp = at ?/> = in ry-models (as r — )• oo, tp^ oc 
and pr^ oc oc tp'^). The function d'^p^/dtp'^ is most conveniently written using the 

intermediate parameter u = 

d^Pn _ T] u^[12 + A(A-r])u + 2(3 -r])u^] 



dtp^ 47r (1 + uf-^ 

^. = ^[i-(i + -)^-1, V^l. ^^^^ 

T] — L 

= ln(l -|- u), r] = 1. 

It is easy to see from equations (17) and (18) that all models with r] < 3 have d'^p^/dtp'^ > 
and hence /,,(e) > 0, so the phase-space density is positive as required. 



For r] > 1 stars are present for all energies between and Cmax = '^/{v ~ 
< r] < 1 stars are present for all e > 0. The DF diverges as e — )• Cmax or e — )• oo 
respectively; we have 



?7 = 1, 



47,5/2 ' 



3 

— ^(emax - e)~\ r] = 3. 



47,3 

(19) 

At small energies, 

M^) ^ as e ^ 0. (20) 

Analytic expressions for the DF are available for a number of ry-models: 



/3W = j^i2(2.)'/^f^+31„g 



1 - (2e)i/2 



.l + (2e)i/2 

f2{e) = ^,/,^3(i_,)5/2 e^/' - [e{l - e)]^/^(3 + 2e - 2A^ + 16e^)} , 

/3/2(^) = 20/.,3(2 - ,)9/2 { 1(3 + 32e - 8e^) sm-^(e/2)^/^ 

^^^^ ~ ^ -(63 + 693e - 5670e^ + 7410e^ - 4488e^ + 1448e^ - 240e^ + 16e^)|, 



28 

m = {-e^erf(.^/^) + ^ erf [(2.)^/^] - ^F(e^/^) + ^F[(2.)^/^]} 



M^) = 29/27,3(2 + ,)9/2 {§(3 - 32e - 8e^)sinh-^(e/2)^/^ 

- ^^^^ ^ ^ (63 - 693e - 5670e^ - 7410e^ - 4488e^ - 1448e^ - 240e^ - 16e^)V 
28 J 

(21) 

where erf(a;) is the error function and F(x) = exp(— x^) exp(t^)c/t is Dawson's integral. 
The expressions for r] = 1 and 2 are due to Jaffe (1983)^ and Hernquist (1990). Analytic 



expressions are also available for other ry-models (e.g. 1] = |, |), but they are even more 
complicated. 



In JafFe's equation (7) the terms involving F[( — E)^^'^] and F[( — 2E)^^'^] should be multiplied by 



-I. 
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2. Phase-space distribution function (mass per unit phase-space volume) /(e) as a function of 
the energy per unit mass — e, for the same values of rj shown in Figure 1. The curves for 77 > 1 
diverge at emax = ^/{v ~ 1)? the depth of the central potential well; for f] < 1 the depth of the 
central potential well is infinite. 



Figure 2 shows the DF /(e) for several ry-models. The DF increases with e for all re- 
models, a result which can be verified analytically and which implies that ry-models are 
stable (BT). 
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Surface brightness Since we assume that the mass-to-hght ratio T is constant, the 
surface brightness at projected radius R is 

m-'^i:^/^. (22) 

Using equation (3) and the substitution u = sin""' (i?/r), the surface brightness of an 
rj-model may be written 



(R + sinuY+i 

Analytic expressions for It]{R) are available for models with integer values of rj: 

1 l-(2- R^ ]X(R] 

^ -3 + (2 + R^)X(R) 

27r(l-i?2)2 ' ' 

_ 2 + 13R^ - {12R^ + 3R'^)X{R) 
~ 47r(l -i?2)3 ' 



(23) 



(24) 



where 



1 1 

COS 



(25) 



^ ' R>1. 



(i?2- 1)1/2 

The expressions for = 1,2 are due to Jaffe (1983) and Hernquist (1990). The surface 
brightness can also be written analytically for half-integer using elliptic integrals; how- 
ever, the expressions are sufficiently complicated that they are not very useful. 

The surface brightness profiles satisfy a recurrence relation: 

Iv+iiR) = - — ^ [R'-'W] ■ (26) 
The asymptotic behavior of the surface brightness is 

T/.(i?) - ^ as i? ^ ^. (27) 



At small radii 



[ln(2/iJ) - |] , „ = 2, (28) 

TT 

1 

2 < ?7 < 3. 



7r(r7 - 2)(r7 - 1) 
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All models with 2 < rj < 3 have a core with constant surface brightness even though 
the volume density diverges as r^~^: the divergence is weak enough that the emissivity 
integrated along the line of sight is dominated by radii of order unity rather than by the 
central cusp. 

For any value of rj the surface brightness may be determined numerically. The numer- 
ical evaluation of the integral in equation (23) is straightforward for i? > 1. For R <^ 1 
numerical integrations are slow to converge because the integrand changes behavior near 
u ^ R (the integrand is u^~^ for R <^ u <^ 1 but / R for u <^ R). A strategy that 
improves the convergence is to integrate numerically the difference between the integral 
in (23) and J^^'^ u^{R + u)~^~'^du^ which has the same behavior for u <^ 1 but can be 
integrated analytically. A subroutine that computes It]{R) is available from the authors. 

The surface brightness profiles of ry-models are plotted in Figure 3(a). Figure 3(b) 
compares the surface brightness profiles to de Vaucouleurs' (1948) empirical law, which 
states that log / is a linear function of i?"'^^. 



Line-of-sight dispersion The velocity dispersion along the line of sight at projected 
radius R is denoted (Tp{R) and is given by 

2 ^ p{r)vl{r)rdr 
' TI{R) (r2 - i?2)i/2 

2 t°° rdr f°° 

= YmJR W^W^l ^''^ 

2 f°° p{r)M{r) 



TI(R) 



Jr f 



where the second line follows from the hydrostatic equilibrium equation (6). Using equa 
tions (3) and (4) we have 

2 Y,iR) 
T/^(i?) 

where 



2,7-3 f 

Jo 



-dx 

,--/2 -3 2 , (31) 

27r Jo (R + sinuy+^'i' 



The integral for Yjj(R) can be evaluated analytically for integer or half-integer rj (Hernquist 
1990 gives the analytic expression for rj = 2). However, the integrals are cumbersome. A 
simpler approach is to use the recursion relation 

Yr,+,/2{R) = [R'-''Y,{R)] , (32) 
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3. Surface brightness of 77-models as a function of radius i?, normalized so that the total luminos- 
ity is unity. The same values of rj are plotted as in previous figures, (a) Log-log coordinates; (b) 
Log of the surface brightness plotted against i?"'^^; in these coordinates de Vaucouleurs' surface 
brightness profile is a straight line. 

together with the analytic expression 

^V.W=|ii-i-3^ + ^5^ + ^^^A-,7J), (33) 

where X is defined in equation (25). 

The asymptotic behavior of the line-of-sight dispersion is given by 

"^,.^^^^15^ asi?^oo. (34) 



12 



The behavior at smaU radii is most easily given in terms of the variable Yj^: 



2 ' 



87ri/2r(3 - riY 

.i^[31n(2/i?)-^], , = |, (35) 



27r(2r] - l)(2r] - 2)(2r] - 3) 



<ri <3. 



From these expressions the behavior of c^(-R) can be deduced using equations (28) and 
(30). For < ?7 < 1, cr^ oc R^~^ diverges; for 1 < < 2 it approaches zero at the center 
(as R^~^ for ?7 < | and as R^~^ for r] > |), and for > 2 it is asymptotically constant at 
the center^. 

The line-of-sight dispersion (Tp{R) is plotted in Figure 4 for several values of rj. The 
subroutine used for these calculations is available from the authors. 



Aperture dispersion The aperture dispersion (Ta{R) is the root-mean-square veloc- 
ity measured through a circular aperture of radius R centered on the galaxy. Clearly 

I{R)al{R)RdR 



(7„ 



I(R)R dR 

The aperture dispersion is only defined for models with > |; otherwise the dispersion 
measured through any central aperture is infinite. A convenient form for numerical evalu- 
ation can be derived from equation (22) and the third line of equations (29): 

J p{r)M{r)r dr - p{r)M{r){r^ - R^f/^dr/r^ 



(7 



3 p[r)r'^dr — p{r)r{r'^ — R^yl'^dr 



i{R) 



^-W-Att p{r)M{r){r^ - R^ fl^dr/ 
^ M -At: p{r)r{r^ - R^y/^dr 



(37) 



where M and W are the total mass and potential energy (eq. 13). For ry-models with ?7 > | 
we have 

_ „2»?-5/ 2 _ ^2^3/2 



a 



6(2ry-l) ^ ,oo r^-\r'-Ry/\ , ' ^^^^ 

At large radii (^^^^(i?) l/[6{2r] - 1)]. The beh avior at small radii is easily deduced from 
equations (28), (30), (35), and (36). For < < |, cr^ oc i?*?"^ as i? ^ 0; for | < < 2, 
cr^ oc R^~^; for > 2, cr^ is asymptotically constant at the center. 

The aperture dispersion is plotted in Figure 5. 

* The leading term of the expression given by Hernquist (1990) for cr|(i?) as i? — > is incorrect. The 
denominator [81n(2/s) — 12]7r should be 41n(2/s) — 6. 
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4. Line-of-sight dispersion function of radius R for the same models shown in earlier 

figures. 

3 Models with a central black hole 

We now consider ry-models that contain a non-luminous central object ("black hole") of 
mass /i. The density p(r) (eq. 3) and surface brightness I{R) (eq. 23) are unchanged, and 
the total mass in stars is still unity. The sole effect of the black hole is to modify the 
potential $(r) (eq. 5) to 

$*(r) = $(r) - ^ (39) 



14 




R 



5. Aperture dispersion Ga (the dispersion measured through an aperture of radius R centered 
on the galaxy) for models with r] = 0.6 (dash-dot line), 1, 2, 3 (solid lines), |-, |- (dashed lines). 
The aperture dispersion is infinite for models with f] < ^. 



(we denote properties of models with a black hole by the superscript *). We restrict 
ourselves to models with 1] < ^ since we show below that models in which the DF depends 
only on energy are unphysical if > | . 
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Velocity dispersion If the velocity-dispersion tensor is isotropic, the radial velocity 
dispersion is given by 



9 * 



(r) = v^ir) + ^ vlir) + ^s\r), (40) 



where u^(r) is given by equations (7) and (8) and 



s 



1 4 6 4 1 



^(r) = r'-^{l + rf^^ - + - + 

— 4 T] — 6 T] — A T] — L T] 

\l+rf ^4:{l + rf 6r{l + rf ^4:r^{l + rf r^(l + r) 



1] — 4 1] — 3 1] — 2 1] — 1 1] 

The cases = 1,2 are special: 

sl(r) = |r"^ - I + |r + 6r^ + 4r^ - 4r^(l + rf ln(l + 1/r), 
sl(r) = \r~^ - I - llr - 15r^ - 6r^ + 6r(l + rf ln(l + 1/r). 

As r ^ 0, 



(41) 



(42) 



which represents the usual r divergence in the mean-square velocity expected near a 



black hole. As r — )• oo, s^{r) — )• |r for all i]. 



Distribution function The DF can be determined analytically when the energy e is 
large (i.e. close to the black hole). We begin by evaluating the integral in equation (16), 

^ r dp dr ^ r^'^ dp du 

Jo dtp* {e-tP*y/^ Jo du[e-tP*{u)Y/^' ^ ' 

where we have changed the integration variable to u = and ^(e) is defined implicitly 
by ?/'*[«(£)] = e. We have 

v-; = + ^ = ^ [1 - (1 + u)'-^] +pu, ^ + 1, 

' r ?7 - 1 ' (45) 

= ln(l -|- u) -|- /iu, ?7 = 1, 



and 



dp^ r] 



3 



When the energy e is sufficiently large, the integral (44) is dominated by the contribution 
from the largest values of u, so that tp* — )• pu, dp^/du — )• [ri(3 — ri)/4:7r]u'^~^ , and 

rj{3-rj) r'^ y'-^du _ rjT{4 - rj) 
^ 47r io (e-puy/^ 47ri/2^3-^r(| - ry) ' ^ ^ 
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Thus 

Since T{x) < for — 1 < x < 0, models with > | are unphysicaL There is a simple 
physical explanation for this constraint: consider a DF consisting of stars at a single energy 
-Eo, f{E) = S(E — -Eo), in the potential of a point mass /i, $(r) = —fJ,/r. The density is 
p{r) = 47r/o°°(5(|t;2 - /x/r - Eo)v'^dv = 2^/^Tr{Eo + ji/rf/^ At smaU radii p{r) r'^/^ 
Thus, if the DF is isotropic, stars of a given energy always produce a density cusp around 
a black hole that rises as r"-'/^ at radii r <^ fj,/\E(,\; ry-models with ?7 > | are not allowed 
because their density cusps are shallower than this limiting value. 

Figure 6 shows the DF for ry-models with several values of the black hole mass /i. 

The stability properties of isotropic spherical stellar systems containing a central black 
hole are not presently known. 



Line-of-sight dispersion The mean-square velocity along the line of sight may be 
written 

^ nMiy, ,49, 

where 



W = 2 / 
Je 



^ , (50) 

= ^R^-' / , ^ 'I dx. 
27r ii (1 + Rxy+^ 

Comparing to equation (31) we find the convenient identity 

Z^(i?) = 2F^/2(i?), (51) 

so that the program used to compute can be used to compute Z^, and analytic expres- 
sions for Z\ and Z2 can be found from equations (32) and (33). At small radii, 

(R) ^7 77^ 1 < ?7 < 2, 



TTjl 



4i?[ln(2/i?) 



31' 
2i 



'7 = 2, (52) 



7,1/2 ( IV 2)r(| - i?7)/i 



8r(3-|ry) 

Figure 7 shows the line-of-sight dispersion profile for several ry-models with black holes 
of varying masses. 
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6. Phase-space distribution function /* as a function of energy e for models with rj = 1, f, 2, 
and I", containing central black holes of mass /i = (dashed line), 0.1, 0.3, and 1. The total mass 
in stars is 1. 



4 Other "simple" models? 

In a strict mathematical sense the ry-models are not particularly simple, since all of the 
quantities we have investigated — radial velocity dispersion, surface brightness, line-of-sight 
dispersion, aperture dispersion, and DF — can be derived by quadratures starting with (say) 
an arbitrary density profile p{r). In a practical sense the models are not the simplest pos- 
sible ones either: except for special values of rj the quadratures for the surface brightness, 
line-of-sight dispersion, aperture dispersion and DF are most easily done numerically. 
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7. Line-of-sight dispersion function of radius R for models with 77 = 1, |-, 2, and |-, 

containing central black holes of mass /i = (dashed line), 0.1, 0.3, and 1. The total mass in stars 
is 1. 



Nevertheless the ry-models share several "user-friendly" features: the density 
mass M^(r), and potential all have simple analytic forms; the mean-square radial 

velocity v^(r) is analytic; and the density can be expressed as a simple function of the po- 
tential (which is useful in Eddington's formula [16]). We have failed to find generalizations 
of or alternatives to the ry-models of comparable simplicity and realism. 

It may be useful to describe a chain of argument that leads naturally to the ry-models. 
In choosing a simple functional form for a galaxy model, the natural place to start is with 
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the potential $(r), since this can be differentiated to yield the mass and density: 



MW = .^-. pM = j^^. (53) 

These operations are simpler when we use the variable u = 1/r: 

MM = --, P(u)=-^. (54, 

Thus it is natural to look for a simple function such that $ — )• —u as u — )• (which nor- 
malizes the total mass to unity); another useful property is that should be invertible, 
so that u can be expressed as a function of $. A suitable choice is = [1 — (1 + 
which leads to ry-models with b = 1 — rj. Other choices are possible but the ones we have 
examined are all less attractive. 

Note from equation (54) that the asymptotic behavior p oc as r — )• oo (/? oc as 
M — > 0) that is common to all the ry-models is a feature of any functional form that can 
be expanded in a Taylor series, $(m) = —u + X]j^2 ^^j'^'' with a2 ^ 0. If the first non-zero 
coefficient is aj,, the density at large radii falls as 



5 Discussion 

The models that we have described have many limitations: they are spherical, whereas 
most elliptical galaxies and spiral bulges are probably triaxial; the DF is isotropic in velocity 
space, an assumption that has no strong justification; and their common surface brightness 
profile at large radii (/ oc R~^) does not adequately represent the diverse behavior of 
the outer envelopes of real galaxies. Despite these shortcomings, the ry-models provide 
useful illustrations of the kinematic and photometric behavior that can be present in the 
central parts of galaxies. They also offer simple "strawman" models for comparison with 
observations of galaxies that exhibit photometric power-law cusps in their centers. 



5.1 Classification of central structure 

As an illustration, let us imagine that a galaxy is observed to have a surface brightness 
profile I{R) which is oc R at small radii, falling more steeply at radii larger than some 
characteristic radius that we shall call the "break radius" . We also assume that the galaxy 
is spherical and its DF is isotropic in velocity space. We distinguish "true" radii, which 
are measured in three dimensions from the center of the galaxy, from "projected" radii, 
which are the components of the true radius vector that are normal to the line of sight. 
Our results imply that three types of central structure are possible, not just for ry-models 
but for any galaxy with these general features: 

(I) flat core structure (7 = 0, corresponding to 3 > rj > 2): A density cusp as steep as 
p oc may be present in this type, but there is no cusp in the surface brightness. 
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which is dominated by stars whose true radii are of order the break radius. When 
observed with hmited resolution, some galaxies in this class may appear to have small 
but non-zero values of 7 (cf. Fig. 3). If there is no central black hole, the depth of the 
central potential well is finite; the line-of-sight dispersion is asymptotically constant 
as i? — )• and is dominated by stars at true radii comparable to the break radius. If a 
black hole is present, the mean-square line-of-sight velocity grows, at a rate between 
oc and R~^'^ depending on 1]; models with ?7 > | are unphysical in this case. 
Hernquist's model [r] = 2) marks the boundary between types I and II. 

(II) weak cusp structure (0 < 7 < 1, 2 > > 1): In this type = 2 — 7; thus the density 
cusp is /? oc r . The surface brightness is dominated by stars whose true radii are 
comparable to the projected radius. If no black hole is present the central potential 
has finite depth and the mean-square line-of-sight velocity approaches zero near the 
center, cr^ oc R"* for 7 < | and cr^ oc R^ for 7 > |. If a black hole is present the 
mean-square velocity rises as oc R~^ . Jaffe's model (7 = 1) marks the boundary 
between Types II and III. 

(Ill) strong cusp structure (1 < 7 < 2, 1 > > 0): In this type, as in the previous one, 
?7 = 2 — 7 and the density cusp is /? oc r . The central potential well has infinite 
depth, and the mean-square line-of-sight velocity diverges near the center, oc R^ 
with no black hole and oc R~^ if a black hole is present. 



5.2 Measuring mass-to-light ratios 

The traditional method of determining the mass-to-light ratio of spherical stellar sys- 
tems with flat cores is known as core fltting or King's method (King 1966, Richstone and 
Tremaine 1986). The determination is based on the formula for the mass-to-light ratio 



Q^liR = 0) 

27rGI{R = 0)i?hb 



T = k——^ ^T^' (55) 



where G is the gravitational constant (set to unity in earlier sections), i?hb is the half- 
brightness radius deflned by /(i?hb) = |-^(0), and k is nearly unity for a wide variety of 
stellar systems (and hence is set to unity in the usual applications of the method). 

Core fltting clearly is not an adequate approach to measuring the mass-to-light ra- 
tios of most ry-models, or of any galaxy with similar central structure. To illustrate the 
difflculties, we evaluate the constant k in equation (55) as a function of the parameter rj. 
For ?7=3, k = 1.048; k rises to 1.091 at r] = 2.64, then sinks to 0.900 hy r] = 2.32. Thus 
for 2.32 < r] < 3 the core fltting formula (55) with A; = 1 is accurate to better than 10%. 
However, for values of rj outside this limited range the accuracy of the core fltting formula 
plummets: at rj = 2.2 we have k = 0.50, at = 2.1 we have k = 0.05, and for rj < 2 the 
formula fails completely since the central surface brightness diverges. 

Other methods must replace core fltting for models with rj < 2.3: 
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(i) The local method is based on the surface brightness I{R) and hne-of-sight dispersion 
(Tp{R) measured at a single radius; dimensional analysis then yields the formula 

^ - "amkr '''' 

where h{ri,R) is a dimensionless constant that is plotted in Figure 8. When < | 
and R is much smaller than the break radius, h is given by the analytic formula 

^=22.-i,3/2r(2-,)r(|-,)' ^^'^^^ ^''^ 

The local method is appropriate when the radius R is large compared to the spatial 
resolution of the observations. 

(ii) The aperture method is based on the aperture dispersion (Ta{R) and the luminosity 
inside an aperture, S{R) = 27r I(R' )R' dR' : 

where g{ri,R) is a dimensionless constant that is plotted in Figure 8. The aperture 
method fails for < | because the aperture dispersion diverges. When | < < | 
and R is much smaller than the break radius, g is given by the analytic formula 

9 22i-V"/2>,r(2-,)r(|->,)' ^ ^ ■'■ ' ' 

The aperture method is appropriate when the radius R is close to the limiting spatial 
resolution. 



5.3 The formation process 

Our discussion so far has focused exclusively on the range of possible equilibrium models, 
without addressing issues of galaxy formation. Different formation processes may favor 
particular values of rj. Several examples of such constraints are known for the case where a 
central black hole is present. If the black hole grows slowly, arguments based on adiabatic 
mvariance imply 7=1 (Peebles 1972, Young 1980, BT), ahhough in this case the velocity- 
dispersion tensor is somewhat anisotropic. If a steady-state distribution of stars around 
the black hole has been established by two-body relaxation, we expect 7 = | for stars of 
equal mass (Bahcall and Wolf 1976, BT). If the formation process leaves no stars bound 
to the black hole, we expect 7 = 0, ?7 = I (Peebles 1972, BT). 

Programs to compute the surface brightness, distribution function, line-of-sight and 
aperture dispersion of ry-models are available by electronic mail from the authors. ST 
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8. The dimensionless quantities g and h (eqs. 56 and 58) that determine the mass-to-light 
ratio from measurements of the velocity dispersion and surface brightness, plotted as a function 
of radius R for models with r] = 0.6 (short dash-long dash curve), 1, 2, 3 (solid curves), |-, |- 
(unlabelled dashed curves). 
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